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Abstract 


Black widows (BWs) are millisecond pulsars ablating their companion stars. The out-flowing material from the 
companion can block the radio emission of the pulsar, resulting in eclipses. In this paper, we construct a model for 
the radio eclipse by calculating the geometry of the bow shock between the winds of the pulsar and companion, 
where the shock shapes the eclipsing medium but had not been described in detail in previous works. The model is 
further used to explain the variations of the flux density and dispersion measure of three BW pulsars (i.e., PSR 
B1957+20, J2055+3829, and J2051—0827) detected by the Five-hundred-meter Aperture Spherical radio 
Telescope. Consequently, we constrained the parameters of the three BW systems such as the inclination angles 
and true anomalies of the observer as well as the mass-loss rates and wind velocity of the companion stars. With 
the help of these constraints, it is expected that magnetic fields of companion stars and even masses of pulsars 
could further be determined as some extra observation can be achieved in the future. 
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1. Introduction 


Millisecond pulsars (MSPs) are rapidly rotating neutron 
stars, with a large number of them being discovered in binary 
systems. It is widely believed that MSPs are spun up via 
angular momentum transfer from their companion star through 
accretion processes (Bhattacharya 1996). Redbacks (RBs) and 
black widows (BWs), widely known as spider pulsars, are sub- 
populations of MSPs that can evaporate their donors 
persistently by their powerful winds. Masses of RB companion 
stars are usually within the range of ~(0.2-0.4)Mz, while BW 
companions with a range of ~(0.02—0.05)M. (Roberts 2011; 
Polzin et al. 2019). The optical observations of RB and BW 
binaries confirmed that the surface of their companion stars 
could be strongly irradiated by the energetic pulsar winds 
(Khechinashvili et al. 2000). Furthermore, observations at 
higher energies also discovered the X-ray emission arising 
from the interaction of the pulsar wind with the evaporation 
material (Kluzniak et al. 1988; Phinney 1988). 

According to observations, the eclipse regions of spider 
pulsar binaries are typically larger than the Roche lobes of the 
companion stars, and therefore the eclipse is not solely caused 
by the companion star’s obscuration, but probably due to the 
absorption of the radio pulse signal by the surrounding plasma 


around the companion star, especially during the beginning and 
ending stages of the eclipse. Since the low-mass companion 
stars are surrounded by a diffuse evaporation halo, the radio 
emission from the MSPs can in principle be eclipsed by the 
evaporation halo as it appears on the line of sight (LOS) of 
observers near inferior conjunctions (Guillemot et al. 2019). 
Such eclipses can happen as long as the orbital inclination of 
the binary is high enough. The first discovered eclipsing BW in 
radio observation is PSR B1957+20 (Fruchter et al. 1988). 
About 10% of its orbital phase was found to be eclipsed, which 
cannot be covered by the range of the Roche lobe of its 
companion and thus the ablated material must replenish the 
eclipse medium to be beyond the Roche lobe (Ray & Loeb 
2017; Polzin et al. 2019). Therefore, it is expected that the 
study of the MSP eclipses can help to probe the property of the 
evaporation material and, furthermore, of the companion stars. 
Meanwhile, a constraint on the orbital parameters of the 
binaries could also be obtained. 

The physical processes of eclipsing pulsars had been investigated 
thoroughly by Thompson et al. (1994), by confronting the 
observational results of PSR B1957+20 and PSR J1748—2446A. 
The correlation between the eclipse duration and the frequency of 
PSR J2215+5135 was studied by Broderick et al. (2016), who 
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suggested that the primary effect of the eclipsing medium is to 
absorb the radio emission rather than scattering it. Polzin et al. 
(2019) reported the radio observations of PSR J2051—0827 in the 
frequency range of 110—4032 MHz and found that the scattering 
and/or cyclotron absorption provides the most promising eclipse 
mechanism. Kudale et al. (2020) further suggested the absorption is 
primarily due to the cyclotron-synchrotron process, by explaining 
the excess dispersion, scattering, and absorption of PSR J1227 
—4853. In any case, it is worth mentioning that cyclotron absorption 
requires a strong magnetic field in the eclipse medium. 

In the last decade, many eclipsing pulsar binaries have been 
discovered, which allows us to study the eclipse mechanisms in 
more detail (Bhattacharyya et al. 2013; Guillemot et al. 2019; 
Nieder et al. 2020). In particular, the Five-hundred-meter 
Aperture Spherical radio Telescope (FAST) has implemented 
plenty of radio observations for several BW pulsars such as 
PSR B1957+20, PSR J2055+3829, and PSR J2051—0827, 
including long-term measurements of their flux density and 
dispersion measure (DM; Wang et al. 2023). Furthermore, 
rotation measure (RM) variation was found in PSR J2051 
—0827. Therefore, at the current stage, it is necessary to 
combine these observational results with the eclipse processes 
to constrain the physical properties of binary systems. 

In the next section, we introduce an eclipse model by taking 
into account the shock interaction between the winds of the pulsar 
and companion star as well as the consequent geometry of the 
shock. Then, we apply the model to fit the light curves and DM 
variations of the three FAST BWs in Section 3. The implications 
of the parameter constraints for the binary systems are discussed 
in Section 4. Finally, a conclusion is given in Section 5. 


2. The Model 


The geometry of the eclipse medium is determined by the 
shock surface between the pulsar wind and stellar wind as well 
as a magnetosphere (Romani & Sanchez 2016; Wadiasingh 
et al. 2017). Recently, Chen et al. (2021) investigated the shock 
interaction of a pulsar with an O/B star in some high-mass 
gamma-ray binaries and estimated the radio emission window 
of these pulsars. The primary difference between the spider 
pulsar binaries from the gamma-ray binaries is that, in the latter 
case, the intense wind from the O/B star can nearly enclose the 
pulsar, whereas in the former case, the stellar wind can only 
blow a conical region. The geometry of the bow shock due to 
the interaction of the spider pulsar wind with the companion 
wind is depicted in Figure 1, which is determined by the 
mechanical balance between the two shocked wind regions. 
The most crucial parameter is the momentum rate ratio between 
the winds as: 


_ Lya/C 


= , 1 
Wea (1) 
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where Lsa is the spin-down power of the pulsar, c is the speed 
of light, Mc is the mass loss rate of the companion star, and vw 
is the terminal speed of this wind. 

Following Chen et al. (2021), when the LOS of observers 
can intersect with the bow shock, the radio emission from the 
pulsar would be suppressed as F,(v) = F,,9-e ™”, where the 
optical depth of the companion wind on the LOS is considered 
to be dominated by the free—free absorption (FFA) process as: 


rv) = J av; T, ne, m)dl, (2) 
p,obs 

where a(v; T, ne, ni) is the FFA coefficient, which is dependent 
on the temperature (T) of the medium and the density of the 
electrons (ne) and ions (nj) , and the value of /, ps is determined 
by the shape of the bow shock. Here, the FFA process rather 
than the cyclotron/synchrotron processes is invoked, because 
the latter processes usually require a strong magnetic field. The 
strength of the magnetic field could in principle be measured by 
detecting the variation of the RM of the radio emission near the 
eclipse boundaries (Polzin et al. 2019; Kudale et al. 2020). For 
the spider pulsar PSR B1957+20, Li et al. (2019) measured its 
RM values and found no evidence for large-scale magnetic 
fields over the egress plasma lensing region. Besides the radio 
suppression effect, the companion wind can further provide an 
extra contribution to the DM of the pulsar radio emission, 
which leads to a DM variation as 


ADM = i: ndl, GB) 
lp, obs 

where ne is the electron number density, which is determined 
by the hydrogen abundance of the wind and the distance to the 
center of the companion star. Specifically, we can relate the 
electron number density with the ion density by ne = nifui/ He» 
where the typical values of the mean ion molecular and electron 
weight are taken as i~ 1.29 and pe~ 1.18, respectively 
(Zdziarski et al. 2010). 

In order to complete the above integrals, we write the 
number density of the companion wind (i.e., the ion density 
n,=Nny) as a function of the radius r to the companion as 
(Waters et al. 1988): 


-2 
r 
Ny (1) = mo 2} > (4) 
Fy 
where the base density nyo at the stellar surface r, can be 
defined as: 
Mc 
Nw 0 = 2 > (5) 
AT Vw himp 

where m, is the mass of protons, vw is the wind velocity at 
surface of companion. Here the mass-loss of the companion is 
assumed to be isotropic, which is viable for the region not very 
close to the companion. In addition, it should be mentioned that 
when the orbital inclination is high enough (i.e., io > 5 — 0), 
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Figure 1. Left: Illustration for the general case of the orbital inclination (i ) and the true anomaly of the LOS (¢,) in spider pulsar binary systems. Right: Description of 


bow shock structure in spider binary systems of the companion star. 


the companion star can appear on the LOS and obstruct the 
radio emission completely, where 0, is the opening angle of the 
companion star with respect to the centroid of the binary. 

On the other hand, the temperature of the stellar wind 
gradually decreases through adiabatic cooling with increasing 
radial distance from the surface of the companion star. As a 
result, the temperature distribution can be empirically 
expressed by a power law as (Kochanek 1993): 


—8 
Lins (=) = Ar-$, (6) 


Fy 


where T, is the effective temperature of the star at r, and 
A= Tr. The value of the index 8 depends on the adiabatic 
index of the wind gas, which could range within ~(2/3-4/3). 
In the following calculations, we simply adopt @ = 2/3. 
Strictly speaking, the stellar surface temperature is actually 
anisotropic due to the irradiation by the pulsar and thus the 
temperature T, involved here can only be treated as an effective 
one. Nevertheless, it should still be noticed that the temperature 
directly involved in the calculation is actually the temperature 
of the wind material far away from the stellar surface, where 
the anisotropy of the temperature is very likely to be not as 
significant as on the stellar surface. 

By confronting the above model with the radio observations 
by spider pulsars, we can constrain the model parameters of the 
binary orbit (the inclination angle i, and the true anomaly of 
observer @,) and of the companion star (Mc, vw, and A). 


3. Fitting the Flux and DM Variations 


3.1. PSR B1957+20 


PSR B1957-+20 is the first MSP discovered in a binary with a 
low-mass companion star, and it is classified to be a BW system. 
The period and eccentricity of the binary orbit are Por» ~ 9.2 hr and 
e<4.0 x 10 > (Arzoumanian et al. 1994). PSR B1957+20 has a 


spin period of 1.607 ms and a spin-down luminosity of 
La~ 1.6 x 10% erg s7! (Fruchter et al. 1988). During the eclipse 
phase, the flux density of PSR B1957-+20 also decreases to nearly 
zero, Which makes it unable to obtain the value of DM. In the left 
panel of Figure 2, we present the observational data centered at 
v~ 1.25 GHz on 2022 September 2, which covers the whole 
eclipsing period and lasts 14,340 s, where the flux density is 
normalized by the pre- and post-eclipse emission (S. Q. Wang 
et al. 2023, in preparation). Solid lines represent fittings of the data, 
and the goodness of fitting is evaluated with the MCMC method. 
The observational constraints on the model parameters are shown 
in the right panel of Figure 2 and the corresponding parameter 
values with lo errors are listed in Table 1. 

The obtained inclination angle of ~85° 06+°:38 indicates the 
orbit of PSR B1957+20 is nearly viewed edge-on, which is 
consistent with the result (i, ~ 85°) given by Johnson et al. (2014) 
through their MSP light curves simulation model. In comparison, 
Reynolds et al. (2007) derived a relatively smaller inclination of 
i, = 65° + 2° from the photometric data of the William Herschel 
Telescope and Hubble Space Telescope. Alternatively, Kandel 
et al. (2021) constrained the inclination to be i, ~ 75°8 + 5°9 by 
using the X-ray emission arising from the stellar wind interactions. 

The mass-loss rate of the companion of PSR B1957+20 is 
constrained to Mc ~ 107115 Mo yr~!, which is also well 
consistent with the result i.e., Mc ~ 107!? Mz yr!) of Polzin 
et al. (2020) who analyzed the 149 MHz observations with an 
inclination of 65° (Reynolds et al. 2007). It is considered that 
the mass loss of the companion star is caused by the 
evaporation due to the irradiation from the MSP and, thus, its 
rate can be estimated by 


2 
Mc = | = (7) 


2 
2v Aorb 


where aor» is the orbital separation of the binary and f is the 
evaporation efficiency. Then, by using Mc ~ 107128 Mo yr}, 
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Figure 2. Left: Normalized flux density and DM variations of PSR B1957+20 during the eclipse phase. The solid lines give the model fittings for the observational 
data. Right: Corner plot showing the posterior distributions of the model parameters. 


Table 1 
Parameters and Fitting Results of PSR B1957+20, PSR J2051—0827, PSR J2055+3829 


Parameters PSR B1957+20 PSR J2051—0827 PSR J2055+3829 
orbital period, Poy, (hr) 9.2 (a) 2.38 (d) 3.1 (g) 
spin-down luminosity, Lg (erg s7!) 1.6 x 10°% (b) 5.5 x 10%” (e) 4.3 x 10”? (g) 
eccentricity, e a075) 4.0 (a) 6.2 (f) 1.4 (g) 
projected binary separation, do SIn io (Ro) 2.7 (c) 1.02 (f) 1.2 (g) 
inclination angle of observer, i, (°) 85.067).38 59.504039 46.804143 
true anomaly of observer, ġo (°) 94.7774 0-05 81,220.02 100.48+9:03 
mass-loss rate, logọoMce (Mo yr’) —12.154985 —11.92*0 01 —11.644 0.54 
wind velocity, vw (10° cm s~') 2.057608 0.5470 34 1.0372" 
temperature coefficient, A (10'° K em?!) 0.234351 2.781 5:88 0.927931 
evaporation efficiency f (107°) 0.235 3.33 4.34 


Reference. (a) Fruchter et al. (1988); (b) Arzoumanian et al. (1994); (c) Lin et al. (2023); (d) Stappers et al. (1996a); (e) Shaifullah et al. (2016); (f) Lazaridis et al. 
(2011); (g) Guillemot et al. (2019). 


we can obtain f= 2.35 x 107° for PSR B1957+20. Further- 
more, we can in principle explore a clue from the mass-loss rate 
to determine the evolutionary stage of the companion star 
(Chen et al. 2013). 


3.2. PSR J2051—0827 


PSR J2051—0827 is the second eclipsing MSP discovered 
after PSR B1957+-20. The MSP is in a 2.38 hr orbit around the 


companion star with an eccentricity of 3 x 1074. The pulsar 
has a spin period of 4.508 ms and a spin-down luminosity of 
6 x 10*ergs~'. Radio observations at low frequency 
(v < 0.6 GHz) show eclipse features covering 10% of the 
orbital period (Stappers 1996). The FAST observations of PSR 
J205 1—0827 centered at v ~ 1.26 GHz was carried out on 2022 
January 14 with a duration of 8940s, covering the entire binary 
orbit (Wang et al. 2023). The observational data are displayed 
in the left panel of Figure 3, which shows the radio emission of 
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Figure 3. The same to Figure 2 but for PSR J2051—0827. 


the pulsar was not completely obscured during the eclipse. 
Furthermore, in comparison with the previous source, the data 
of PSR B1957+20 has a much larger scatter, which makes it 
impossible to obtain a good fitting with a simple model. A 
possible explanation for the large fluctuation of the data is that 
the stellar wind in this source is very unstable and in- 
homogeneous. 

In any case, by according to a rough fitting to the data of 
PSR J2051—0827, we can still get a constraint on the 


inclination angle as 59° 50+0;49, In the previous work of 


Dhillon et al. (2022), this angle was found to be 55° 9+4;8 by 
fitting the symmetrical HiPERCAM light curves with a direct- 
heating model of the pulsar to the companion star. Meanwhile, 
Stappers et al. (2001) estimated it to be ~40° by using the 
asymmetry optical light curves of the companion, where the 
gravitational distortion of the companion and the irradiation of 
the pulsar wind are taken into account. The mass-loss rate of 
the companion Mc ~ 107!!9? Mo yr~! is also found to be 
close to the previous estimate by Polzin et al. (2019), which 
determines an evaporation efficiency of 3.33 x 107°. 


3.3. PSR J2055+3829 


PSR J2055+3829 was first discovered by the SPAN512 
survey conducted with the Nangay Radio Telescope, which has 
a spin period of 2.08 ms and a spin-down luminosity of 
La~ 4.3 x 10° ergs”! (Guillemot et al. 2019). As another 
BW pulsar, PSR J2055+3829 is in a binary with a very low 


mass companion of a tight orbit. S. Q. Wang et al. (2023, in 
preparation) conducted observations on PSR J2055+3829 for 
11,640 s centered at 1.25 GHz on 2021 August 30, and found a 
clear eclipse with a duration of 1467+ 12s, which is about 
13% of the orbital period. The observational data are shown in 
the left panel of Figure 4, in comparison with the modeling 
curves. 

The fitting results indicate a relatively moderate inclination 
angle of i, ~ 46° 8041:33, which is between the edge-on orbit 
(i ~ 90°) assumed by Guillemot et al. (2019) and the inclination 
of 26° estimated from the mass function with an upper limit on 
the companion mass. The mass loss rate and the velocity of the 
companion wind are constrained to be Mc ~ 107!!- Mo yr! 
and 1.03 x 10°cms_!, respectively, indicating an evaporation 
efficiency of f= 4.34 x 107°. 


4. Discussion and Implications 
4.1. RM Measurements and Magnetic Fields 


For PSR J2051—0827, Polzin et al. (2019) utilized the 
linearly and circularly polarized flux to investigate the possible 
presence of magnetic fields through the measurements of RM, 
Faraday delay, and depolarization. As a result, the average 
magnetic field parallel to the LOS was constrained to be 
within the range of 10 “*G< By <S 10° G. Very recently, 
Wang et al. (2023) detected a regular decrease of the RM of 
PSR J2051—0827 in the egress of the eclipse. i.e., the RM 
gradually change back to normal when the LOS moves away 
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Figure 5. A comparison of the model-predicted RM variation of PSR J2051—0827 during the eclipse phase with the observational data (solid circles), where a typical 
magnetic field structure is adopted as presented in Equation (9). Different values for the parameter B, and the inclination of the magnetic axis @m are taken as labeled 
in the legend, while the other model parameters are taken as listed in Table 1. The observational data are taken from Wang et al. (2023). 


from the eclipse medium, as shown by the data points in 
Figure 5. This measurement provided direct evidence for the 
existence of magnetic fields in the eclipse medium. 

In theory, based on the obtained constraints on the orbital 
and companion wind properties, we can in principle predict the 
possible variation of RM of the spider pulsars by 


3 lo @) 
ARM = —“_ f neBy dl, (8) 
2rmé ct Ip.obs 


by adopting a specific description of the magnetic field in the 
stellar wind (Bw). Typically, we can express the magnetic field 
as a combination of the radial and toroidal components as 


(Bosch-Ramon & Khangulyan 2011): 


By(r, o) = Bye; + Byes (9) 

with 
B, ~ B,(r,/T), (10) 
By ~ B, (Vw, /Yw) (r/r), (11) 


where B, is the magnetic field strength at r, and the surface 
rotation velocity of the companion star is assumed to be 
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vw, ~ 0.lvẹ(r,/r). Finally, it needs to be pointed out that the 
symmetric axis of the above magnetic field could be not 
parallel to the angular momentum of the binary, but probably 
has an inclination angle of ¢,, relative to the normal of the 
orbital plane, which would also influence the RM variation. 
Here the binary is assumed to be tidally locked. 

By using the above formulae, we plot the theoretical RM 
variation curves in Figure 5 for different magnetic field 
strengths and different magnetic inclinations, where the 
orbital and wind parameters are taken as listed in Table 1 
for PSR J2051—0827. However, by comparing the theoretical 
curves with the observational data, it is unfortunately found 
that the model always cannot explain the data. An obvious 
difference is that the data is delayed by a time of about 
0.015% of the orbital period relative to the model prediction. 
This indicates that it is necessary to consider a more 
complicated structure for the magnetic field, in particular, 
around the bow shock (e.g., Phinney 1988). Due to the fact 
that the magnetic field near the bow shock can be significantly 
compressed by the shock and, furthermore, the distribution of 
the magnetic field may not completely trace the bow shock. In 
particular, the LOS could be scanned by the tail of the bow 
shock for a long period because the shock tail can be deflected 
due to the orbital motion. In this tail region, the magnetic field 
could be very high so that a significantly high RM 
contribution can be expected, whereas the corresponding 
radio absorption and DM contribution are negligible. 
Additionally, the orbital motion of the companion would not 
only have an influence on the shock tail but would also lead 
the symmetric axis of the bow shock to deflect from the 
connecting line of the binary. The deflection angle of the bow 
shock can in principle be estimated by € ~ arctan(vc/vw), 
where vc is the orbital velocity of the companion (Zabalza 
et al. 2013). For PSR J2051—0827, we can get €~47°5, 
which is much larger than the other two BW binaries. 
Therefore, the consideration of these complications about the 
magnetic field configuration around the bow shock is 
definitely needed in future modeling works as more RM data 
can be obtained. 


4.2. The Companion Surface Temperature and Pulsar 
Masses 


As the radius of the companion star of PSR J2051—0827 had 
been measured to be r, = 0.139t89!1 R> by Dhillon et al. 
(2022), we can roughly derive its surface temperature to be 
T, = 6093+134 K from the obtained parameter value of 
A = 2.7818% x 10!° K cm2/3, This temperature is higher than 
the one given by Stappers et al. (1996b) as T, = 4000-4700 K, 
which was estimated from the observed R — J color according 
to the cool M-type dwarf relations of Bessell. For PSR J2055 
+3829, if we use the Roche lobe radius Rg ~ 0.14 Rẹ to 
approximate the companion radius, which is reasonable for an 
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Figure 6. The relationship between the masses of PSR B1957+20, PSR J2055 
+3829, and PSR J2051—0827 and their companions. 


appropriate mass ratio between the pulsar and the companion 
(Guillemot et al. 2019), we can estimate the surface 
temperature of the companion to be 2016*35K from the 
parameter value A = 2.7808 x 10!° K cm?/3, However, for 
PSR B1957+20, the combination of the obtained parameter 
A = 0.23193! x 10!°K cm?/3 with the possible radius 
r, œ 0.25Rẹ of the companion (van Kerkwijk et al. 2011) 
would lead to a temperature as low as T, ~ 343+!3 K, which is 
much lower than that found by Reynolds et al. (2007) and is 
probably unrealistic. This indicates that the value of 8 = 2/3 
may not be always appropriate for the companions or 
Equation (6) for a constant index is not strictly correct, 
especially for the wind very close to the stellar surface. In 
addition, we would like to emphasize again that the 
temperature of the companion should in fact be anisotropic 
due to the irradiation of pulsar wind and the temperature 
obtained from the fitting of the radio flux and DM is only an 
effective value, which might be regarded as an average 
temperature. 

In any case, the determination of the companion surface 
temperature would in principle help us to calculate their 
emission luminosity and, furthermore, to estimate the masses of 
the companions by according to an empirical mass—luminosity 
relationship. Meanwhile, the determination of the inclination 
angle of the BW binaries enables us to fix the mass function of 
the binaries as 


(12) 


On ? (ap sin ig)? _ (Mcsin io)? 
Rb G (Mp + M.) ` 


where M, and Mc are the masses of the pulsar and companion, 
respectively. To be specific, the Mp — Mc relationships of the 
three FAST BW binaries are displayed in Figure 6. Then, if the 
companion masses can indeed be obtained from the 
temperature measurement, it is expected that the pulsar masses 
would finally be constrained. 
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5. Summary and Conclusion 


The study of pulsar eclipse processes occurring in BW 
systems can provide a more accurate perception of binary orbits 
and their companions. It is proposed that the radio emission 
from the pulsar is absorbed by the outflows of its companion. 
The interaction of the pulsar wind with the stellar wind forms a 
bow shock, and the stellar outflows in the shock-wrapped 
companion side cause the radio eclipse when the star moves 
around inferior conjunctions. 

We establish an eclipse model of wind interaction scenario 
and fit the radio flux density and DM variations of three BW 
pulsars detected by FAST, namely PSR B1957+-20, PSR J2055 
+3829, and PSR J2051—0827. We constrain some of the 
orbital parameters as well as the wind parameters of companion 
stars under the eclipse model. For PSR B1957+20, we found 
that the on has a relatively high orbital inclination 
(i, = 85° (oe 38) and a true anomaly of the observer of 
dy = 94° T1002 The mass loss rate and the wind velocity of 
the low-mass companion are Mc ~ 107!>!5 Mo yr! and 
vw ~ 2.05 x 108cms~'. While for PSR J2051—0827, we 
obtain the orbital parameters such as its orbital inclination of 

= 59? 5040:40 and the true anomaly of the observer as 


i = 81227. ape the companion’s mass loss rate and wind 


velocity are constrained to be Mc ~ 107!!9? Mo yr! and 
vy ~ 0.54 x 108 cms, respectively. Similarly, for PSR J2055 
+3829, we use the eclipse model to constrain the orbital 


inclination as i, = 46° 80tl; ey and the true anomaly of the 


observer as ġ, = 100° 4870s ees and the mass loss rate and the 
wind velocity of the companion are constrained to be 
Mc ~ 1071164 Mo yr™! and vw ~ 1.03 x 108cms™', respec- 
tively. The above fitting results indicate that our eclipse model 
can be used to constrain the orbital parameters as well as the 
wind parameters of the companions in BW binary systems. 

In addition, our eclipse model still needs further refinement, 
as for the evidence of the magnetic field in the RM variation 
detected in the late stage of the eclipse. We suggest that it is 
not dominated by the magnetic field of the companion star but 
by the interaction between the magnetic field of the pulsar 
wind and the bow shock, which relies on a more accurate 
description of the magnetic field around the bow shock. 
Furthermore, if the mass and radius of the companion can be 
given directly by optical observations, we can also make an 
estimate of the mass of the pulsar in the BW systems by 
eclipse. 
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